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ABSTRACT 

We conducted Very Long Baseline Array (VLBA) observations of seven nearby narrow-line 
Seyfert 1 (NLS1) galaxies at 1.7 GHz (A18cm) with milli-arcsccond resolution. This is the first sys- 
tematic very long baseline interferometry (VLBI) study focusing on the central parsec-scale regions 
of radio-quiet NLSls. Five of the seven were detected at a brightness temperature of > 5 x 10 6 K and 
contain radio cores with high brightness temperatures of > 6 x 10 K, indicating a nonthermal process 
driven by jet-producing central engines as is observed in radio-loud NLSls and other active galactic 
nucleus (AGN) classes. VLBA images of MRK 1239, MRK 705, and MRK 766 exhibit parsec-scale 
jets with clear linear structures. A large portion of the radio power comes from diffuse emission com- 
ponents that are distributed within the nuclear regions (< 300 pc), which is a common characteristic 
throughout the observed NLSls. Jet kinetic powers limited by the Eddington limit may be insufficient 
to allow the jets escape to kiloparsec scales for these radio-quiet NLSls with low-mass black holes of 
< 10 7 Mq. 

Subject headings: galaxies: active — galaxies: jets — galaxies: Seyfert — radio continuum: galaxies 
- galaxies: individual (MRK 110, MRK 705, MRK 1239, NGC 4051, MRK 766 
(NGC 4253), MRK 783, and AKN 564) 



1. INTRODUCTION 

Narrow-line Seyfert 1 galaxies (NLSls) are identified 
by their optical properties: a flux ratio [O III]/H/3<3, 
and Balmer lines only slightly broader than forbid- 
den lines (Osterbrock & Pogge 1985), defined as 
FWHM(H/3)<2000 km s^ 1 (Goodrich 1989). NLSls 
make up ~ 20% of an optically selected active galac- 
tic nucleus (AGN) sample from the Sloan Digital Sky 
Survey (SDSS; e.g., Zhou et al. 2006). Recently, 7-ray 
emissions have been detected from several very radio- 
loud NLSls by the Large Area Telescope onboard the 
Fermi Gamma-ray Space Telescope (Abdo et al. 2009a, b; 
Foschini 2011); it has been argued that these NLSls may 
represent a third class of 7-ray emitting AGNs along with 
blazars and radio galaxies. The detected 7-rays strongly 
suggest the presence of blazar-like phenomena as a con- 
sequence of relativistic beaming effects due to a fast jet 
aligned close to our line of sight (Blandford & Konigl 
1979). Moreover, the presence of highly relativistic jets 
in radio-loud NLSls has been implied on the basis of ra- 
dio flux variability observations (Zhou et al. 2003; Doi 
et al. 2006, 2011; Yuan et al. 2008), very-long baseline 
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interferometry (VLBI) imaging (Doi et al. 2006, 2011; 
Giroletti & Panessa 2009; D'Ammando et al. 2012), and 
modeling of the spectral energy distributions (Zhou et al. 
2007; Yuan et al. 2008; Abdo et al. 2009b). Kiloparsec 
(kpc)-scale radio morphology is associated with several 
radio-loud NLSls, including two 7-ray emitting NLSls 
(Anton et al. 2008; Gliozzi et al. 2010; Doi et al. 2012). 
However, these powerful jet phenomena are observed in 
only a small part of the NLS1 populations. In fact, the 
fraction of radio-loud objects 7 in the NLS1 population is 
significantly low (Zhou & Wang 2002; Stepanian et al. 
2003; Greene et al. 2006; Zhou et al. 2006; Yuan et al. 
2008), ~ 7% with R > 10 and - 2.5% with R > 100 (Ko- 
mossa et al. 2006) , compared with normal Seyfert galax- 
ies and quasars (~ 10%— 15%; e.g., Ivezic et al. 2002). 
Similarly, kpc-scale radio emissions are less common in 
NLSls than in broad-line AGNs (Doi et al. 2012). The 
radio jet characteristics of a huge amount of the radio- 
quiet NLS1 population behind these extremes are little 
known. 

NLSls as a class are considered to have weak jet activi- 
ties. A radio-imaging study at arcsecond resolution using 
the Very Large Array (VLA) revealed that nearby (radio- 
quiet) NLSls, in most cases, are unresolved radio sources 
[< 300 parsec (pc)] with modest radio luminosities (10 20 ~ 
10 23 W Hz -1 ) and steep radio spectra (a < —0.5, where 
S v oc v a , and S u is the flux density at the frequency 
v). Thus, NLSls are not remarkably different from nor- 
mal Seyfert galaxies at radio wavelengths ( Ulvestad et al. 
1995, see also Moran 2000; Zuther et al. 2011), although 
they exhibit many distinctive optical-X-ray properties. 

7 The radio loudness R is defined as the ratio of the 5 GHz 
radio flux density and optical B-band flux density. A threshold of 
_R = 10 separating radio-loud and radio-quiet objects (Kellermann 
et al. 1989) is frequently used. 
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In addition to the unusually narrow B aimer lines, NLSls 
exhibit several extreme properties such as strong permit- 
ted Fe II emission lines (Boroson & Green 1992), rapid X- 
ray variability (Pounds et al. 1995; Leighly 1999b), steep 
soft X-ray spectra (Wang et al. 1996; Boiler et al. 1996; 
Leighly 1999a), and frequently observed blue-shifted line 
profiles (Zamanov et al. 2002; Boroson 2005; Leighly & 
Moore 2004) . It has been suggested that these properties 
are related to high-mass accretion rates close to the Ed- 
dington limit (Boroson & Green 1992; Brandt & Boiler 
1998; Sulentic et al. 2000; Mineshige et al. 2000) on rela- 
tively low-mass black holes (~ 10 5 -10 7 - 5 M Q ; Peterson et 
al. 2000; Hayashida 2000; Grupe & Mathur 2004; Zhou et 
al. 2006). A high- accretion rate on a low- mass black hole 
supports the radio-quiescence of the NLS1 class, accord- 
ing to many studies using large AGN samples: a positive 
correlation between radio loudness and the black hole 
mass (Laor 2000; Lacy et al. 2001; Dunlop et al. 2003; 
McLure & Jarvis 2004; Metcalf & Magliocchctti 2006) 
and an anti-correlation between radio loudness and the 
accretion rate (Ho 2002; Lacy et al. 2001; Maccarone et 
al. 2003; Greene et al. 2006). Similarly, Merloni et al. 
(2003) discovered "a fundamental plane of black hole ac- 
tivity" in three-dimensional space consisting of the radio 
luminosity, X-ray luminosity, and black hole mass (rang- 
ing from 10 1 to 10 10 M Q ) for samples of X-ray binaries, 
our Galactic center, low-luminosity AGNs, Seyfert galax- 
ies, and quasars. These tendencies in jet activity can be 
attributed to nonlinear dependences on the black hole 
mass and accretion rate (Heinz & Sunyaev 2003). Thus, 
the weak radio emissions of NLSls as a class can be un- 
derstood in the framework established for other AGN 
classes. 

However, the picture of a nonthermal jet origin for the 
observed radio emissions of NLSls has not been fully es- 
tablished yet, particularly for the radio-quiet subclass, 
because of the difficulty of VLBI observations of weak 
radio sources. Previously, a nonthermal jet origin for 
radio-quiet quasars was established by a series of VLBI 
detections (Blundell & Beasley 1998; Ulvestad et al. 
2005). The first VLBI study focusing on the NLS1 class 
at milli-arcsecond (mas) resolutions examined the very 
radio-loud narrow-line quasar PMN J0948+0022 using 
the Very Long Baseline Array (VLB A) (Doi et al. 2006). 
A Doppler-beaming effect of highly relativistic jets is re- 
quired to explain the observed extremely high brightness 
and rapid flux variability in the VLBA images. Sub- 
sequent 7-ray detection by Fermi/LAT followed (Abdo 
et al. 2009a). Doi et al. (2007) detected five additional 
radio-loud NLSls using the Japanese VLBI Network and 
established the picture of a nonthermal process in the 
radio-loud NLS1 subclass. Gu & Chen (2010) reported 
compact radio structures for three of the five radio-loud 
NLSls using VLBA imaging. Doi et al. (2011) found ev- 
idence of pc-scale jet structures in all the five radio-loud 
NLSls and pointed out (1) the presence of relativistic 
nonthermal jets, which are essentially the same as those 
of other radio-loud AGN classes and (2) the existence 
of both apparently radio-loud NLSls owing to Dopplcr 
beaming and intrinsically radio-loud NLSls. Several 7- 
ray emitting NLSls have been observed in survey pro- 
grams of VLBI radio sources: SBS 0846+513 (Kovalev 
et al. 2007; D'Ammando et al. 2012), 1H 0323+342 (Zhou 
et al. 2007), and PKS 1502+036 (Fey & Chariot 2000, 



see also Dallacasa et al. 1998 for a study of this object 
as a compact radio galaxy candidate) and are now being 
monitored by VLBA in the MOJAVE 8 program. On the 
other hand, there have been only a handful of images 
of radio-quiet NLSls: MRK 766, AKN 564 (Lai et al. 
2004), NGC 5506 (Middelberg et al. 2004), NGC 4051 
(Giroletti & Panessa 2009), and MRK 1239 (Orienti 
& Prieto 2010) as part of a study of Seyfert galaxies 
(not focused on NLSls). A systematic VLBI study of 
radio-quiet NLSls is crucial for understanding the nature 
of their radio-quiescent central engines with the lowest- 
mass black holes at the highest accretion rates in AGN 
classes. Furthermore, radio-loud NLSls should be dis- 
cussed as extreme examples on the basis of more solid 
knowledge of radio-quiet objects that are the majority in 
the NLS1 population. 

In the present paper, we report VLBI imaging of seven 
nearby radio-quiet NLSls. In Section 2, we describe the 
sample selection. In Section 3, we describe the observa- 
tions and imaging processes. The observational results 
are presented in Section 4. We discuss the observed radio 
properties of the radio-quiet NLSls in Section 5. Finally, 
we summarize our studies in Section 6. The observed 
properties of individual radio-quiet NLSls are described 
in Appendix A. Throughout this paper, a ACDM cos- 
mology with H a = 71 km s _1 Mpc" 1 , fl M = 0.27, and 
Oa = 0.73 is adopted (Komatsu et al. 2009). 

2. SAMPLE 

A nearby radio-quiet NLS1 sample consisting of seven 
objects was made for this VLBI imaging study (Table 1). 
We selected from the 205 NLSls listed in "Quasars and 
Active Galactic Nuclei," 10th Ed. (Veron-Cetty & Veron 
2001) and the 64 NLSls in "Spectrophotometric Atlas of 
NLSls" (Veron-Cetty et al. 2001). We processed the fol- 
lowing selection criteria step-by-step: (1) redshift z < 0.1 
(129 sources survived), (2) inclusion in the VLA Faint 
Images of the Radio Sky at Twenty-Centimeters (FIRST; 
Becker et al. 1995), 03 April version (60 sources sur- 
vived), (3) existence of a radio counterpart with an offset 
of <1" in the FIRST image (33 sources survived), (4) a 
peak intensity of more than 5 mJy beam -1 (13 sources 
survived), and (5) availability of a radio position mea- 
sured at a resolution of < 0'.'5 from literature or archival 
data (to enable processing by VLBA correlator as de- 
scribed in detail in Section 3.1; 7 sources survived). One 
of the selected NLSls, MRK 1388, has sometimes been 
classified as a type 1.9 or 2 Seyfert (e.g., Goodrich 1989) 
with an unresolvable broad-line component (Osterbrock 
1985; Veron et al. 1997); hence, we excluded it from 
our sample. In addition, we included AKN 564, one 
of the best-studied NLSls, although it does not prop- 
erly satisfy criterion (2); we substituted VLA archival 
data (AT 149) obtained at <~ 5 arcsec resolution as equiv- 
alent to the FIRST data. As a result, seven sources were 
selected. They are genuine NLSls according to the op- 
tical definitions and also other aspects such as the Fe II 
line strength, soft X-ray spectrum, and black hole mass 
[Columns (3)-(9) in Table 1]. All seven NLSls are ra- 
dio quiet by the definition of radio loudness using the 
nuclear optical and total radio flux densities [i?*; Col- 

8 Monitoring of jets in AGNs with VLBA experiments (Lister et 
al. 2009) 
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TABLE 1 

Nearby radio-quiet NLS1 sample for VLBA Observations. 



Name 
(1) 


z 

(2) 


FWHM(H/3) [O 
(km s- 1 ) 
(3) 


iii]/H/3 K4570 Ref. 
(4) (5) (6) 


r 

(7) 


log M BH 

(M e ) 
(8) 


Ref. 
(9) 


Qtotal 
°1.4GHz 
(mjy) 

(10) 


(ii) 




(PC) 

(12) 


Ref. 
(13) 


MRK 110 


0.0353 


1760 


0.71 


0.26 


1 


2.59 


6.75 


4 


10.1 


1.8 


2100 


8 


MRK 705 


0.0292 


1790 


0.41 


0.36 


2 


2.53 


6.72 


5 


8.5 C 


2.1 


< 55 


9 


MRK 1239 


0.0199 


1075 


1.29 


0.63 


2 


4.06 


5.89 


5 


62.8 


0.69 


46 


10 


NGC 4051 


0.0023 


1120 


0.55 


0.25 


2 


3.04 


6.13 


4 


14.7 d 


0.62 


330 


11 


MRK 766 


0.0129 


1630 


0.35 


1.83 


2 


2.63 


6.80 


6 


39.7 


0.63 


65 


12 


MRK 783 


0.0672 


1655 


2.29 


0.49 


2,3 


1.30 a 


7.16 


7 


33.2 


2.0 


< 380 


13 


AKN 564 


0.0247 


865 


0.92 


0.67 


2 


3.4 b 


6.75 


4 


29.1 


0.44 


320 


9 



Note. — Col. (1) source name; Col. (2) redshift; Col. (3) FWHM of H/3 line width of broad-line component; Col. (4) flux ratio of 
[O III]A5007 to H/3; Col. (5) flux ratio of Fe II multiplets in the range 4434-4684A to H/3; Col. (6) reference for optical-line emissions: 
Col. (7) ROSAT soft X-ray photon index at 0.1-2.4 keV (Pfefferkorn et al. 2001); Col. (8) black hole mass; Col. (9) reference for the 
black hole mass; Col. (10) NRAO VLA Sky Survey (NVSS; ~ 45") flux density at 1.4 GHz as total flux density; Col. (11) radio loudness 
defined as the ratio of fc-corrected 5 GHz radio to optical B-band flux densities (Kellermann et al. 1989). We calculated the radio loudness 
using the NVSS flux density and SDSS 7DR PSF magnitude (Abazajian et al. 2009) extinction-corrected by assuming a Balmcr decrement 
of H«/H/3 = 3.1 for all sources, except AKN 564 for which we used a V-band nuclear magnitude (Mullaney & Ward 2008). We assume 
a = -0.7 for radio (except for MRK 1239, MRK 766, and AKN 564; see Appendix A) and a = -0.61 for optical (Grupe et al. 2010) for 
the k corrections. Col. (12) maximum extent of radio structure from nucleus in VLA images at arcsecond resolution; Col. (13) reference 
for the radio size. 

References. — (1) Grupe et al. 2004, (2) Veron-Cetty et al. 2001, (3) Kovacevic et al. 2010, (4) Kaspi et al. 2000, (5) Ryan et al. 2007, 
(6) Done & Gierlihski 2005, (7) Wang & Lu 2001, (8) Kukula et al. 1998, (9) Schmitt et al. 2001, (10) this paper (Appendix A), (11) Ho 
& Ulvestad 2001, (12) Ulvestad et al. 1995, (13) Ulvcstad & Wilson 1984. 

a A large uncertainty (±0.72) 

b Bollcr ct al. 1996 

c Flux density from FIRST (~ 5") at 1.4 GHz to avoid the con- 
tamination from an unrelated source ~ 30" south (Gallo ct al. 
2005) 

d Flux density from VLA A-configuration (~ 1'.'2) at 1.4 GHz (Ho 
& Ulvestad 2001), because the NVSS emission is contaminated by 
the contribution of a host galaxy (Baum ct al. 1993) 



umn (11) in Table 1]. They all, except for MRK 110, are 
compact radio sources covering < 300 pc in VLA images 
[Column (12) in Table 1]. 

3. OBSERVATIONS AND DATA REDUCTION 

3.1. Sky Positions 

The ficld-of-vicw (FOV) in this study is limited for sta- 
tistical reasons. Searching for emission with a detection 
limit of 4cr requires a FOV limited to ~ 32,000 points, 
corresponding to a radius of ~ 100 times the synthesized 
beam on the sky. The expected size of the synthesized 
beam is about 11 masx5 mas aligned nearly north-south 
at the declinations of our sample: the FOV 9 is supposed 
to be ~ 17 1x075. Hence, we attempted to determine the 
positions for correlation (the phase-tracking centers) to 
sub-arcsec accuracy during sample selection (Section 2). 

For each NLS1, we obtained the radio position with an 
uncertainty of ~ O'.'l in images from the VLA A- array 
configuration in the archive (AM384, AK394, AB973, 
and AL410) and/or in the literature (Ulvestad & Wilson 
1984; Ulvestad et al. 1995; Kukula et al. 1998; Nagar et 
al. 1999; Thean et al. 2000; Schmitt et al. 2001). We used 
the National Radio Astronomy Observatory (NRAO) 
SCHED software when it was necessary to convert a B1950 
VLA position to J2000. The phase-tracking centers used 
in our VLBI observations are listed in Table 2. 

3.2. Observations 

9 The bandwidth smearing for a nominal continuum mode and 
the time-average smearing for 2-sec integration in the VLBA corre- 
lator restrict the FOV to ~ 7" and ~ 15", respectively (Thompson 
ct al. 1986; Wrobel 1995). 



We observed the seven NLSls on May 09, 2005 at 
1.7 GHz using ten antennas of the VLBA at the NRAO 
(observation code BD106). A left circular polarization 
was received at a center frequency of 1.667 GHz with 
a total bandwidth of 32 MHz. We observed in phase- 
referencing mode, which allows us to derive calibra- 
tion parameters for instrumental and atmospheric effects 
from the observations of nearby strong compact radio 
sources ("calibrators;" Beasley & Conway 1995). The 
phase-reference cycle periods were 5 min, and the total 
on-source time was about 40 min for each target. 

3.3. Data Reduction 

Standard calibration procedures for VLBA phase ref- 
erencing were applied during data reduction using the 
Astronomical Image Processing System (MPS; Greisen 
2003). Inappropriate earth-orientation parameters were 
corrected using the AIPS task CLC0R. Amplitude calibra- 
tion using a-priori gain values together with system tem- 
peratures measured during the observations was applied; 
the accuracy is typically 5%. The ionospheric dispersive 
delay was corrected using a global ionospheric model and 
the task TECDR. The time and spectral domains were care- 
fully inspected and data flagging was performed; there 
was no significant radio interference in the data. 

Conservatively, we first searched for emission compo- 
nents with peak intensities of > 5er of the noise on im- 
ages with natural weighting in a FOV of l'/l x 0'.'5, 
and also greater than 6cr in a FOV of 370 x 173 (Sec- 
tion 3.1). At this stage, we detected three sources: 
MRK 110, MRK 1239, and MRK 783. They initially 
appeared at positions offset by less than several tens of 
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TABLE 2 

Parameters of phase-reference VLBA observations. 



Name 



(1) 



Phase-tracking 
R.A. ( h m s ) 
(2) 



center (J2000) Rcf. 
Decl. (° ' ") 

(3) (4) 



Astromctric position (J2000.0) 
R.A. ( h m s ) Dccl. (° ' ") 



(5) 



(6) 



Reference 
Source 
(7) 



A9 

(°) 
(8) 



MRK 110 
MRK 705 
MRK 1239 
NGC 4051 
MRK 766 
MRK 783 
AKN 564 



09 25 12.8500 
09 26 03.2763 
09 52 19.0996 
12 03 09.6105 

12 18 26.5186 

13 02 58.8430 
22 42 39.3397 



+52 17 10.388 
+12 44 03.659 
-01 36 43.553 
+44 31 52.680 
+29 48 46.471 
+16 24 27.748 
+29 43 30.915 



09 25 12.8476 +52 17 10.387 



09 52 19.0980 

b 



13 02 58.8438 

b 



-01 36 43.546 



+ 16 24 27.777 



J0932+5306 1.40 
J0931+1414 1.94 
J0948+0022 2.16 
J1155+4555 1.98 
J1217+3007 0.33 
J1311+1658 2.09 
J2236+2828 1.85 



Note. — Col. (1) source name; Col. (2)-(3) position of phase-tracking center used in VLBA correlator; Col. (4) reference used to 
determine the phase tracking center; Col. (5)-(6) astrometric positions determined by our VLBA observations; Col. (7) reference calibrator; 
Col. (8) separation angle between target and reference calibrator. 

References. — (1) average of the results of Ulvestad & Wilson 1984, Kukula et al. 1998, and VLA archival data AL410; (2) AB973; 
(3) average of the results of Thean et al. 2000, Ulvestad et al. 1995, AK394, and AM384; (4) Hagiwara et al. 2003; (5) average of the results 
of Kukula et al. 1995 and Thean et al. 2001; (6) Ulvestad & Wilson 1984. 



a Position is not available because of detection in self-calibrated 
image 
b Non detection 

mas from the phase-tracking centers. After correcting 
the map centers to the emission peaks, we deconvolved 
the images using the task IMAGR (CLEAN) and mea- 
sured the astrometric positions of the emissions using 
the task JMFIT (Table 2). The images of negative detec- 
tion sources MRK 705, MRK 766, and AKN 564 showed 
quadrangular patterns resembling beam patterns. There 
were several emission candidates with intensities of Aa- 
5cr in the images. The patterns seemed to contribute to 
increasing the image noise to levels significantly higher 
than the expected thermal noises. This suggested that 
significant emission does exist, but its power could be 
scattered owing to imperfect calibration of the phase- 
referenced array. The quality of the phase referencing 
is usually more or less degraded by the uncertainty in 
the correlator models, mainly regarding the residual at- 
mospheric (troposphere and ionosphere) delay for each 
antenna's site. Hence, we attempted self-calibration in 
phase with a signal-to- noise ratio threshold of 2.5 using 
the task CALIB. Astrometric information was lost in the 
self calibration procedure. We defined the success of self- 
calibration according to whether the solutions were found 
as a smooth-phase gradient for all the antennas through 
more than half the observing time on a target. The self- 
calibration was successful and emissions were detected 
for MRK 705, MRK 1239, MRK 766, and MRK 783. We 
smoothed and interpolated the solutions. The deconvo- 
lution and self-calibration algorithms were interactively 
applied a few times. As a result, NGC 4051 and AKN 564 
were not detected with a detection limit of 5 times the 
image noise. 

We used the Difmap software (Shepherd 1997) to 
make final images using the CLEAN algorithm from 
the calibrated visibilities (additional self-calibration was 
not performed in Difmap). We used step-by-step uni- 
form, natural, and (it, t>)-tapered weighting functions 
for CLEAN to retrieve both compact and diffuse com- 
ponents as much as possible. Consequently, we ob- 
tained self-calibrated images for four targets (MRK 705, 
MRK 1239, MRK 766, and MRK 783) and a direct image 
for one target (MRK 110). 



4. RESULTS 

Five of the seven radio-quiet NLSls were detected in 
the VLBA observations. The results of the image anal- 
yses are shown in Table 3. The images with natural 
weighting for the detected sources are shown in Figure 1. 
We also show images convolved with a restored beam 
size of 10 pc in Figure 2. We describe the details for 
individual radio-quiet NLSls in Appendix A. 

4.1. Brightness Temperatures 

The images of the five detected radio-quiet NLSls show 
a compact component (core) with high brightness tem- 
peratures of Tb > 6 x 10 7 K (Table 3), which were cal- 
culated from 



T B = 1.8 x 10 9 (1 + z)- 



S v 



(1) 



in Kelvin at the rest frame, where z is the redshift, 
is the flux density in mJy at the observing frequency v 



in GHz, and 



and 



in mas are the fitted full 



widths at half-maximum of the major and minor axes of 
the component, respectively (cf. Ulvestad et al. 2005). If 
a source was unresolved, we adopted half the beam sizes 
as the upper limits to the source sizes (</> < 0/2). 

The typical image noise of ~ 0.41 mJy beam -1 cor- 
responds to T B ~ 1.7 x 10 6 K. Thus, VLBI imaging 
works as a high-brightness filter. Even diffuse emissions 
of <~ 3a (the lowest contour level in the images) imply 
Tb ~ 5 x 10 6 K, which is too high to attribute to any 
stellar origin. We conclude that all the detected VLBA 
emissions from these radio-quiet NLSls are relevant to 
the activity of AGNs. 

4.2. Source Structures 

The five detected sources have compact emissions with 
high brightness temperatures. Furthermore, we can rec- 
ognize significant structures in the images of MRK 1239 
(Figure 1), MRK 705, and MRK 766 (Figure 2). A lin- 
ear structure that extending to cither side of a brightness 
peak was clearly seen in MRK 1239. The visible radio 
extent was ^100 mas (corresponding to ~ 40 pc), which 
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Fig. 1. — VLBA images with natural weighting for detected nearby radio-quiet NLSls at 1.7 GHz in the region of 150 X 150 mas 2 . 
Contour levels are separated by factors of v2 beginning at 3 times the rms noise (Table 3). Negative and positive contours are represented 
by dashed and solid curves, respectively. Half-power beam sizes are given in the lower left corners. Angular scales corresponding to a linear 
scale of 10 pc are given in the lower right corners. 



was consistent with a slight elongation found in a VLA 
archival image at sub-arcsec resolution [Figure Al (b) in 
Appendix A] . MRK 705 showed a one-sided linear struc- 
ture ~ 45 mas (~ 26 pc) east of the core. No larger 
structure has been found in VLA images (Appendix A). 
MRK 766 showed a signature of a northern elongation 
of ~ 50 mas (~ 13 pc) in the restored image (Figure 2). 
The VLBA images for the other two sources (MRK 110 
and MRK 783) exhibit apparently compact emission. 

The VLBA-detected flux density was significantly 
lower than the VLA total flux density for all the sources, 



indicating low fractions of high-brightness components 
(> 5 x 10 6 K) relative to lower-brightness components 
(Section 4.1). In this paper, the radio core fraction, 
/i tqhzi is defined as the ratio of the flux density from 
the central 10 pc (in diameter) to the total flux density 
at 1.7 GHz in the rest frame. At the distances to our 
samples, 10 pc typically corresponds to ~ 20 mas [Col- 
umn (10) in Table 3]. Our VLBA baselines can retrieve 
> 95% of the emissions as a correlated flux density at the 
minimum (projected) length of ~ 210 km for the 20-mas 
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Fig. 2.— 

radio size. Hence, the fluxes from the central 10 pc were 
practically measurable [Column (11) in Table 3]. On the 
other hand, we cannot rule out the possibility that emis- 
sions with much larger angular scales were partially or 
scarcely detected ("resolved out"). Hence, we adopted 
the total flux densities from the NRAO VLA Sky Sur- 
vey (NVSS; Condon et al. 1998) at 1.4 GHz as listed 
in Table 1; the contribution of extended component was 
calculated by subtracting the < 10 pc contribution at 
1.4 GHz from the NVSS value at the rest frame, and 
then, converted to the value at 1.7 GHz, assuming spec- 
tral indices of compact and extended components (see, 
Appendix A and Table 3 in detail). Consequently, the 



core fractions were typically less than 10%, up to 24% at 
most [Column (12) in Table 3]. Because the entire radio 
structures are < 300 pc in the VLA images [Column (12) 
in Table 1], a significant portion of the emissions should 
originate in the regions of ~ 10-300 pc in the form of 
low-brightness components detected (> 5 x 10 6 K) and 
not detected (< 5 x 10 6 K) in the VLBA images. In 
fact, we recognized a clear enhancement in the visibility 
amplitude of MRK 1239 at baselines of < 3 MA [Fig- 
ure Al (a) in Appendix A] and marginal enhancements 
also in MRK 766 and AKN 564 (Appendix A). These 
are evidences for the presence of a significant amount of 
emission components with sizes of > 70 mas, correspond- 
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TABLE 3 
Results of VLBA observations. 



Name rVLBIcorc oVLBIcorc oVLBItotal n \ 7j~. p~A tZ flTT, l?< iu P c W 

IName J 1.7GHz ^l.rGHz ^l.TGHz CT "maj frnin B »10pc ^i.7GHz Jl. 

(mjy beam -1 ) (mjy) (mJy) (mjy beam -1 ) (mas) (mas) (°) (xl0 7 K) (mas) (mjy) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


MRK 110 


1.1 ± 0.1 


1.2 ± 0.2 


1.1 


0.128 


10.9 


4.4 


-12.7 


> 6.5 


14.4 


1.1 


0.13 


MRK 705 


1.8 ±0.2 


1.8 ±0.3 


3.3 


0.144 


11.2 


4.9 


-0.8 


> 8.7 


17.3 


1.8 


0.24 


MRK 1239 


2.4 ± 0.2 


4.8 ± 0.5 


20.0 


0.169 


11.3 


4.5 


-2.5 


7.4 


25.1 


13.1 


0.23 


NGC 4051 


< 0.7 


< 0.7 


< 0.7 


0.143 


11.2 


4.3 




< 1.0 


121 


< 0.7 


< 0.05 


MRK 766 


1.3 ±0.1 


1.6 ± 0.2 


2.5 


0.103 


11.2 


4.4 


-4.1 


8.4 


38.4 


2.4 


0.07 


MRK 783 


1.4 ±0.1 


1.3 ±0.2 


1.3 


0.129 


10.8 


4.4 


-6.7 


> 7.7 


7.9 


1.3 


0.04 


AKN 564 


< 1.0 


< 1.0 


< 1.0 


0.208 


11.8 


4.5 




< 1.3 


20.4 


< 1.0 


< 0.04 



Note. — Col. (1) source name; Col. (2) peak intensity of core component; Col. (3) flux density of core component, measured by a 
single Gaussian profile fitting using the AIPS task JMFIT; Col. (4) total flux density retrieved by CLEAN components; Col. (5) image rms 
noise on blank sky; Col. (6)— (8) synthesized beam pattern in major axis, minor axis, and position angle, respectively; Col. (9) brightness 
temperature of core component, calculated using Equation (1); Col. (10) angular scale in mas corresponding to a linear scale of 10 pc; 
Col. (11) flux density from central 10 pc at 1.7 GHz in the observer frame; Col. (12) radio core fraction, defined as the ratio of the flux 
density from the central 10 pc to the total flux density at 1.7 GHz in the rest frame. The extended flux component was estimated by 
subtracting the contribution from the central 10 pc [Col. (11)] from the NVSS total flux density at 1.4 GHz (Table 1), assuming the core's 
spectral index available in the literature or a = —0.5 for values unavailable in the literature and a = —0.7 for the extended component. 
See Appendix A for individual objects. 



ing to > 20-30 pc for these sources. 

5. DISCUSSION 

5.1. Compact Emission Components 

Our VLBA imaging provides the first insight into the 
pc-scale region of several nearby radio-quiet NLSls in a 
radio band. All five detected radio-quiet NLSls have a 
compact radio component with a high-brightness tem- 
peratures of > 6 x 10 7 K. This is an evidence of a non- 
thermal process driven by jet-producing central engines in 
radio-quiet NLSls, as well as radio-quiet Seyfert nuclei 
(e.g., Preuss & Fosbury 1983; Neff & de Bruyn 1983) and 
radio-quiet quasars (Blundell & Beasley 1998; Ulvestad 
et al. 2005). 

It is notable that the central components make only 
minor contributions compared to the VLA total flux den- 
sities (i.e., low values of 7 (j Hz ; Section 4.2). A previous 
study using the VLBA reported that all but 1 of 39 radio- 
selected low-luminosity AGNs were detected, and they 
exhibited core-dominated structures in most cases (Na- 
gar et al. 2005, see also Falcke et al. 2000; Anderson et 
al. 2004). However, the radio-quiet NLSls in our sample 
present a contrasting aspect, although both of the AGN 
subclasses are weak radio sources in the nearby universe. 
The radio spectra are also quite different: flat or mildly 
inverted spectra in low-luminosity AGNs (Nagar et al. 
2001; Anderson et al. 2004; Doi et al. 2005, 2011, 2013) 
but steep spectra in NLSls (Moran 2000; Zhou & Wang 
2002). Such flat or inverted spectra would be accounted 
for by optically thick synchrotron emitters at the base 
of the jets (e.g., Falcke & Biermann 1999). The scale- 
invariant jet model predicts that the radio luminosity 
from the jet bases depends on cx Mg^ 12 in most accre- 
tion scenarios, where Mbh is the black hole mass (Heinz 
& Sunyaev 2003). In fact, the known radio sources of 
low-luminosity AGNs tend to have large black holes of 
M BH ~ 10 8 -10 9 M Q (e.g., Nagar et al. 2002b) in contrast 
to values of ~ 10 6 -10 7 M Q for the radio-quiet NLSls in 
our sample (Table 1) and the majority of nearby NLSls 
(e.g., Doi et al. 2012; Grupe & Mathur 2004; Zhou et al. 
2006). Hence, compared with low- luminosity AGNs, the 



radio-quiet NLSls possess less efficient radio cores. 

Moreover, the radio-quiet NLSls also exhibit signifi- 
cantly lower core contributions than radio-loud NLSls 
(Doi et al. 2011); core enhancement due to the Dopplcr 
beaming effect is a natural explanation in some cases 
(Doi et al. 2007; Gu & Chen 2010). In our VLBA im- 
ages with limited dynamic ranges, we found no clear 
evidence of highly relativistic jets, such as unusually 
high-brightness temperatures requiring Doppler boosting 
(Kellermann & Pauliny-Toth 1969; Readhead 1994), for 
the radio-quiet NLSls. We will present a comprehensive 
study of the core activity of both radio-quiet and radio- 
loud NLSls in another paper (Doi et al. in preparation). 

5.2. Extended Emission Components 

Past VLA surveys of NLSls found a combination of two 
radio properties: compact sources (< 300 pc) and very 
steep spectra (Ulvestad et al. 1995; Moran 2000; Zhou & 
Wang 2002). This had indicated that the radio sources 
should be localized, but physically not so compact that 
they become optically thin in synchrotron radiation. Our 
VLBA observations have confirmed this picture at the pc 
scale. 

As an example, MRK 1239 is only slightly resolved in 
VLA images at arcsecond resolutions and is revealed to 
exhibit very diffuse structures within 40 pc in the VLBA 
image. Less than half of the total fluxes were detected 
in the VLBA image and originate in low-brightness com- 
ponents of ~ 5 x 10 6 K (Appendix A). We examine the 
possibility of the thermal process (free-free emission) for 
radio-emitting plasma of 15.2 mJy (the VLBA flux den- 
sity of the extended component for MRK 1239). The 
brightness temperature of 5 x 10 6 K gives a physical 
temperature of < 5 x 10 6 K. According to the Gaunt 
factors of - 10 and - 0.4 at - 1 GHz (Gayet 1970) and 
~ 1 kcV (Sutherland 1998) for plasma at ~ 10 6 K, the 
expected X-ray luminosity would be > 4 x 10 44 erg s _1 
at 0.1-2.4 keV (Rybicki & Lightman 1979), much higher 
than 1 x 10 42 erg s _1 observed by ROSAT toward 
MRK 1239 (Zimmermann et al. 2001, see a similar dis- 
cussion for radio-quiet quasars by Laor & Behar 2008; 
Stcenbrugge et al. 2011). It is difficult to reconcile a 
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thermal process with the observed low-brightness radio 
emissions. The remaining lower brightness components 
(< 10 6 K) that were below the VLBA sensitivity are 
presumably nonthermal jets as well because the position 
angle of the VLA-resolved structure of MRK 1239 [Fig- 
ure Al (b) in Appendix A] is consistent with the direction 
of the two-sided jet-like structure in the VLBA image. A 
steep spectrum in total flux [Figure Al (c) Appendix A] 
also strongly supports a nonthermal synchrotron origin 
for the diffuse emissions. 

MRK 705 and MRK 766 also exhibit clear morpho- 
logical evidences of jets in the VLBA images. A large 
portion of the total flux density comes from diffuse com- 
ponents (f 1 7 Gjj z = 0.24 and 0.07, respectively), which 
is also the case for MRK 1239 (/ 1 1 % c Hz = 0.23). The 
other NLSls also showed low core fractions (f x ^guz ~ 
0.10). NGC 4051 was not detected in our observation 
with fiTQnz < 0.05 but was detected by Giroletti & 
Panessa (2009) at 1.6 GHz by very sensitive VLBI imag- 
ing. The radio morphology is symmetric and shows two- 
sided lobe- like emissions separated by 20 pc and centered 
at a core of sub-mJy, which is aligned with larger-scale 
emissions in the VLA images (Appendix A). There is 
an additional radio-quiet NLS1 that was previously ob- 
served by VLBI. The radio core of NGC 5506, a hidden 
NLS1 candidate (Nagar et al. 2002a), is surrounded by 
several diffuse components in the pc-scale region; a bridg- 
ing morphology indicates a physical connection between 
the core and the diffuse components via jets (Middelberg 
et al. 2004). Furthermore, a low-surface-brightness halo 
~ 350 pc in diameter enshrouds the central features in a 
VLA image (Orienti & Prieto 2010), and only one-third 
of the single-dish flux was detected by VLBI (Middelberg 
et al. 2004). Thus, a large portion of the radio luminosi- 
ties come from diffuse emission components distributed 
within the nuclear regions (< 300 pc), which is a general 
property of non-thermal jets in radio-quiet NLSls. 

These diffuse components are putative reservoirs of dis- 
sipated jet kinetic powers. For our samples (Mbh ~ 10 6 
10 7 M Q : Table 1) and the majority of nearby NLS1 pop- 
ulations (< 10 7 M Q ; e.g., Zhou et al. 2006; Doi et al. 
2012), the jet kinetic powers limited by the Eddington 
limit on the low-mass black holes may be insufficient to 
cause the jets to escape to kpc scales outside the host 
galaxies. Indeed, only a small number of known NLSls 
exhibit kpc-scale radio morphology (Anton et al. 2008; 
Gliozzi et al. 2010; Doi et al. 2012); all of these NLSls 
possess relatively high- mass black holes of > 10 7 M Q and 
estimated jet kinetic powers of > 10 44 erg s _1 , and show 
signatures of rclativistic jets in most cases (Doi et al. 
2012). On the other hand, even the most luminous radio 
source in our radio-quiet sample (MRK 705) possesses an 
extended radio power of 10 23 5 W Hz -1 at 1.4 GHz, sug- 
gesting an estimated jet kinetic power of 10 43 5 erg s~ 4 
(by using the relation of O'Sullivan et al. 2011). These 
values are below the thresholds of the Fanaroff-Riley 
type I/II (Fanaroff & Riley 1974) dichotomy at a radio 
power of - 10 24 - 6 W Hz" 1 (Zamfir et al. 2008) and a 
kinetic power of ~ 10 44 erg s" 1 (Punsly & Zhang 2011) 



that enables escape in the form of supersonic radio lobes 
(Kawakatu et al. 2009; see the discussion in Doi et al. 
2012 for NLSls with kpc-scale radio structures). Fur- 
thermore, advancing radio lobes would be more likely to 
become disrupted by the dense interstellar medium in the 
nuclear region (e.g., De Young 1993) of spiral hosts, as 
seen in normal Seyfert galaxies (Orienti & Prieto 2010). 
The diffuse components found in our NLS1 sample may 
also be relics of such disrupted structures. On the basis of 
the VLBI-observed radio-quiet and radio-loud NLSls, we 
will present a study of their jet kinetic powers estimated 
from the extended radio components (Doi et al. in prepa- 
ration). 

6. SUMMARY AND CONCLUSIONS 

Our VLBA imaging study at 1.7 GHz and mas reso- 
lution provide the first insight into the pc-scale region 
of several nearby radio-quiet NLSls in a radio band. 
Five of seven sources (MRK 110, MRK 705, MRK 1239, 
MRK 766, and MRK 783) were detected at brightness 
temperatures of > 5 x 10 6 K and contain radio cores with 
high-brightness temperatures of > 6 x 10 7 K, indicating a 
nonthermal process driven by a jet-producing central en- 
gine, in common with other AGN classes. MRK 1239 
exhibits a two-sided jet-like structure ~ 40 pc in ex- 
tent. MRK 705 and MRK 766 exhibit one-sided linear 
jet structures of ~ 26 pc and ~ 13 pc, respectively. On 
the other hand, two sources (NGC 4051 and AKN 564) 
were not detected at the sensitivity of our observations. 

We defined the radio core fraction, f 1 7 ojj z , as the ratio 
of the flux density from the central 10 pc (in diameter) to 
the total flux density at 1.7 GHz in the rest frame. We 
found typical values of f 1 7 ojj z < 0.1 and a maximum 
of 0.24 in our sample. While radio sources are appar- 
ently compact at arcsec resolution (< 300 pc), our VLBI 
study indicated that a large portion of the radio lumi- 
nosities come from diffuse emission components that are 
distributed within the nuclear regions. Their central en- 
gines with low-mass black holes of ~ 10 6 -10 7 M Q may 
supply jets with insufficient kinetic powers to escape to 
the outsides of the host galaxies. 
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cility operated by the University of Manchester at Jodrell 
Bank Observatory on behalf of STFC. We used NASA's 
Astrophysics Data System Abstract Service (ADS) and 
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are operated by the Jet Propulsion Laboratory. This 
work was partially supported by a Grant-in- Aid for Sci- 
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APPENDIX 
INDIVIDUAL RADIO-QUIET NLS1S 
MRK 110 

Our VLBA image at 1.7 GHz (Figure 1) showed a single component of 1.1 mJy, which retrieved 13% of the total 
flux density in the rest frame on the basis of the 1.4-GHz NVSS flux density and a spectral index of a = —0.7 between 
1.4 and 8.4 GHz (Kukula et al. 1998). Because we could not make self-calibration on this source owing to insufficient 
signal-to-noise ratios, we cannot rule out the possibility that the VLBA image was affected by some coherence loss 
in phase referencing. VLA A-array images showed a secondary component 3'.'0 (2.1 kpc) to the north of a central 
component (Kukula et al. 1998, see also Miller ct al. 1993). 

MRK 705 

The self-calibrated VLBA image of MRK 705 at 1.7 GHz showed a compact core of 1.8 mJy, and a possible linear 
structure that extends ~ 45 mas (~ 26 pc) to the east (Figures 1 and 2). An unresolved source of 8.5 mJy was found 
in the FIRST image and accompanies a secondary source of 27.5 mJy, which has an unrelated infrared counterpart 
(Gallo et al. 2005). The VLA A-array image at 8.46 GHz showed an unresolved (< 55 pc) source of 2.0 mJy (Schmitt 
et al. 2001), which is consistent with the extent of the VLBA structure. A single component of 5.7 mJy beam -1 was 
seen at 1.7 GHz in an archival image from the Multi-Element Radio Linked Interferometer Network (MERLIN) at 
150 mas resolution (MERLIN proposal code: Y/32); a large portion of the MERLIN fluxes was not retrieved in the 
VLBA baselines. Because the emissions of MRK 705 were detected after self-calibration (Section 3.3), its astrometric 
position could not be obtained. 

MRK 1239 

MRK 1239 is an unresolved (< 80 pc) source showing a steep spectrum in VLA A-array images (Ulvestad et al. 
1995; Thean et al. 2000). However, our VLA A-array image from a good uw-coverage at 8.5 GHz from the archival 
data AM384 showed a significantly elongated structure of ~ 0'.'12, corresponding to ~ 46 pc, at a position angle of 
47° (Figure Al (b)). 

The self-calibrated VLBA image of MRK 1239 at 1.7 GHz shows a low-brightness (2.4 mJy beam -1 or less) linear 
structure of ~ 100 mas (~ 40 pc) consisting of several components (Figure 1), which is consistent with the elongated 
size seen in the VLA image. The visibility amplitudes were clearly enhanced in baselines at < 3 MA [Figure Al (a)]. 
A visibility-based model fitting with a single- Gaussian source structure resulted in a major axis of 62 mas (~24 pc) 
at a position angle of 47°. The flux density of the Gaussian model was 21.5 mJy, consistent with the integrated flux 
density of 20.0 mJy in the VLBA image (Table 3). Orienti & Prieto (2010) also reported a VLBA image of MRK 1239 
at 1.7 GHz that showed two components of 10.2 mJy in total aligned at a position angle similar to our result. The 
contribution from the central 10 pc was 23% of the total flux density at the rest frame on the basis of the 1.4-GHz 
NVSS flux density and assumed spectral indices of a = —0.5 for the core and a = —0.54 (below) for the extended 
components, respectively. 

The presence of diffuse emission distributed within ~ 50 pc can also be inferred from our multi-frequency VLA 
B-array (quasi-simultaneous) observations [AD489; Figure Al (c)]. A spectral index of a = —0.52 was found at 1.4- 
5 GHz; steeper spectra appear at higher frequencies, a = —0.94 at 5-8.5 GHz and a = —1.64 at 8.5-15 GHz (S v oc v a ). 
Such an apparent spectral break also appeared in the VLA A-array (not simultaneous) observations at 1.4-8.4 GHz 
(not simultaneously) by Ulvestad et al. (1995). Moran (2000) pointed out that NLSls tend to have significantly steeper 
radio spectra (a ~ 1.1-1.2) than classical Seyferts on the basis of VLA A-configuration observations at 1.4 and 8.4 GHz. 
They noted that the strong radiation field around highly accreting disks of NLSls can steepen the synchrotron spectra 
owing to inverse-Compton cooling of relativistic electrons in jets. However, our measurements for MRK 1239 from 
uv-tapered images of a common resolution resulted in a simple power-law spectrum with a = —0.54. This provides 
evidence for resolution effects. The significantly steeper spectra generally seen in radio-quiet NLSls might have been 
artificial, which may be related to the dominance of diffuse components (Section 5.2). 

NGC 4051 

NGC 4051 is one of two NLSls not detected in our VLBA observation. To determine the phase-tracking center, we 
reanalyzed the VLA A-array 22-GHz image of Hagiwara et al. (2003), which included nuclear H2O masers confined 
to a region of 0'/04; the VLA observation used the ICRF calibrator J1153+4931 as a reference position. Giroletti & 
Panessa (2009) reported VLBI detection of NGC 4051 using the European VLBI Network (EVN) at both 1.6 and 
5 GHz at a position that differed by only < 10 mas from our phase-tracking center (Table 2). Hence, it was unlikely 
that the nucleus of NGC 4051 was outside of the VLBA's FOV. 

The EVN image at 1.7 GHz showed a core (0.45 mJy) and two symmetric sub-components (0.73 and 0.67 mJy) 
separated by ~ 20 pc with respect to the core with a maximum intensity of ~ 0.34 mJy beam -1 (Giroletti & 
Panessa 2009). Such weak emissions are consistent with our negative detection because they are 3a in the VLBA 
image. A MERLIN image at 1.6 GHz showed a triple structure with an extent of 0'.'8 (38 pc) at a position angle of 
73° (Christopoulou et al. 1997), which corresponds to the structure seen in the EVN image. The VLA A-array image 
at 1.4 GHz showed larger structure extending up to ~ 4" (~ 330 pc) from the nucleus (Ho & Ulvestad 2001; Baum et 
al. 1993); the position angle has changed to 41°, perpendicular to the host-galaxy major axis (135°). The NVSS image 
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shows very diffuse emission associated with the spiral arms and disk at a scale of arcmin, indicating stellar origins 
(Baum et al. 1993). Hence, we used a flux density of 14.7 mJy in the VLA A-array image at 1.4 GHz presented by Ho 
& Ulvestad (2001) as the total flux density from an AGN contribution to derive the core fraction, assuming a = —0.7 
(Christopoulou et al. 1997). 

MRK 766 

The self-calibrated VLB A image with natural weighting at 1.7 GHz showed marginal structures extending to the 
southwest (PA^ —142°) and north (PA^ 0°) from the brightest component (Figure 1). In a restored image, diffuse 
emissions extending up to ~50 mas (corresponding to ~13 pc) north were clearly seen (Figure 2). We recognized an 
enhancement in the visibility amplitude of MRK 766 at baselines of < 5 MA; this suggests the presence of a significant 
amount of diffuse component with the size of > 10 pc. Earlier VLA A-array images showed a slightly resolved structure 
(< 100 pc) at a position angle of ~ 30° (Ulvestad et al. 1995; Nagar et al. 1999; Thean et al. 2000). A MERLIN 
image at 1.7 GHz clearly showed a radio structure extending to 0'.'25 (65 pc) at a position angle of 34° (Thean et al. 
2001), which is consistent with the southwest elongation in our VLBA image with natural weighting. The VLBA flux 
density was 2.5 mJy, which retrieved 7% of the total flux density at the rest frame on the basis of the 1.4-GHz NVSS 
flux density and a = -0.72 between 1.4 and 8.4 GHz (Thean et al. 2001). Lai et al. (2004) reported a VLBI detection 
of MRK 766 at a relatively high intensity 7.1 mJy beam -1 with a synthesized beam of ~2 mas at 5 GHz, which 
was inconsistent with a MERLIN archival image at 5 GHz (MERLIN proposal code: 92NOVA) that showed diffuse 
emission of only 1.8 mJy beam , this discrepancy might be because of variability. Because of the negative detection 
before self-calibration (Section 3.2), the astrometric position could not be obtained from our VLBA observation. 

MRK 783 

Our VLBA image at 1.7 GHz (Figure 1) showed a single component of 1.3 mJy, which retrieved 4% of the total flux 
density at the rest frame on the basis of the 1.4-GHz NVSS flux density and assumed spectral indices of a = —0.5 and 
—0.7 for the VLBA and diffuse components, respectively. The VLA A-array image at 4.9 GHz showed an unresolved 
(< 380 pc) source of 3.5 mJy (Ulvestad & Wilson 1984). 

AKN 564 

AKN 564 is one of two NLSls not detected in our VLBA observation. The VLA A-array image at 8.46 GHz showed 
a triple structure (3.1, 3.0, and 0.9 mJy) extending 320 pc to the north (Schmitt et al. 2001; Moran 2000). A MERLIN 
archival image showed somewhat diffuse emissions of only 1.04 mJy beam -1 at 5 GHz with a 40 mas beam, which 
seemed to be consistent with our negative detection. On the other hand, a mas-scale structure with 3.7 mJy beam -1 
had been detected at 5 GHz by conventional fringe fitting (i.e., not phase referencing) using a much more sensitive 
VLBI array of ^2 mas resolution at 5 GHz by Lai et al. (2004). However they mentioned that they were not confident 
about the structure. THe phase-tracking center of our VLBA observation was determined from the VLA archival data 
AB973 obtained with the A-array at 8.4 GHz (Table 1) using the ICRF position of the phase calibrator J2236+2828 
and differed by only ^260 mas from that of Lai et al. (2004). Hence, it was unlikely that the nucleus of AKN 564 was 
outside of the VLBA's FOV. Some amount of diffuse emissions could be in the VLBA FOV because the correlated 
amplitude was enhanced at baselines of < 5 MA, and the image noise was significantly higher than the expected 
thermal noise level (Section 3.3). Hence, we concluded that AKN 564 is dominated by diffuse components at pc scales, 
which are responsible for our negative detection. The spectral index between 1.4 and 4.9 GHz is a — —0.74 (Ulvestad 
et al. 1981). 




Projected baseline length (MA) RIG1IT ascension (bwsoi Frequency (GHz) 

Fig. Al. — (a) Plot of visibility amplitude of MRK 1239 from VLBA data. Visibilities integrated every 60 min are divided into 30 bins 
over the range of the baseline lengths, (b) VLA images of MRK 1239 with VLA A-array at 8.5 GHz. Contour levels are 3<r X —1, 1, 2, 4, 
8, 16, 32, and 64, where a = 0.16 mJy beam -1 is the image rms noise. The FWHM of the synthesized beam, 0'.'27 X 0'.'22 at PA = 7.6°, is 
illustrated in the lower left corner, (c) Radio spectra of MRK 1239 with VLA B-array. Filled circles connected with solid lines represent 
measurements from VLA images with uniform weighting. Open squares connected with dashed lines represent measurements from tapered 
images with a common beam size using the same data. 
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